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Abstract 



This paper is concerned with the analysis of a class of impacting systems of relevance 
in applications: cam-follower systems. We show that these systems, which can be mod- 
i-G ' elled as discontinuously forced impact oscillators, can exhibit complex behaviour due to 

C^ , the detachment at high rotational speeds between the follower and the cam. We propose 

that the observed phenomena can be explained in terms of a novel type of discontinuity- 
induced bifurcation, termed as corner-impact. We present a complete analysis of this 
bifurcation in the case of non-autonomous impact oscillator and explain the transition 
to chaos observed in a representative cam-follower example. The theoretical findings 



^r\ ' are validated numerically. 

o' 



1 Introduction 



r — ' Recently, much research effort has been spent to analyse the dynamics of piece- 

f^ , wise smooth dynamical systems with impacts lf5l l4T]| . These systems arise in 

many areas of engineering and applied science. A typical example is that of 
mechanical systems characterised by structural components with displacement 
/^ • constraints. Examples include bouncing or hopping robots, systems with back- 

lash or friction, gears, vibro-impacting mechanical devices 0. 

Cam-follower systems are a particularly important class of mechanical sys- 
tems with displacement constraints widely used for the operation of various 
machines and mechanical devices [30]. Usually, their purpose is to actuate 
valves or other mechanisms through the movement of a follower forced by a 
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rotating cam. For example, all types of automated production machines, in- 
cluding screw machines, spring winders and assembly machines, rely heavily 
on this kind of systems for their operation. One of the most common appli- 
cation is to the valve train of internal combustion engines (ICE) |18|, where 
the effectiveness of the ICE is based on the proper working of a cam-f ollo^ver 
system. A schematic of a single valve for a typical pushrod type engine is pre- 
sented in Figure[T] Here, the cam rotation results in a linear motion imparted to 
the valve. The valve spring in the system provides the restoring force necessary 
to maintain contact between the components. 

To guarantee that the follower moves as required, it is important in appli- 
cations to carefully design the cam profile. Different cam geometries are used 
in practice ranging from circular cams to highly complex cam profiles. In gen- 
eral, there is now a large variety of alternative methods to select the cam profile. 
For example, by using constrained optimization algorithm, it is possible to use 
splines to obtain the cam geometry from the desired motion that the cam is 
required to impart on the follower (see for examples [9] and ill6|). This often 
means that while the cam has a continuous displacement profile, it might have 
discontinuities in its acceleration tSTl . 

It has been observed that, as the cam rotational speed increases, the fol- 
lower can detach from the cam. This causes the onset of undesired behaviour 
associated to impacts taking place between the follower and the cam. For ex- 
ample, in automotive engines this phenomenon can seriously deteriorate the 
engine performance as the valves can close with abnormally high velocity and 
even bounce off the seat (valve floating and bouncing) | |2Tll37l[T0l . To avoid 
this phenomenon, a large spring force and preload are applied to the follower 
l34l . This often causes an increase in the contact force, which induces higher 
stresses possibly leading to early surface failure of the parts. The resulting high 
friction valve train reduces the efficiency of the engine system ['39]. 

In general, cam-follower systems can be thought of as impact oscillators 
with moving boundaries |201l5Ul[T5ll40l . While the dynamics of impact oscilla- 
tors with continuous forcing has been the subject of many papers in the exist- 
ing literature (see for example |32| , flTl , ||6l[7|), the possible intricate bifurcation 
behaviour of impact oscillators with discontinuous forcing was discussed only 
recently, as for example in [8]. It was proposed that discontinuously forced 
oscillators can show a novel bifurcation phenomenon unique to their nature 
which was termed as corner-impact bifurcation. Namely, in fS) the d5mam- 
ics are studied of an impact oscillator forced by a discontinuous sinusoidal 
forcing of the form f{t) = A\ sm{cut)\. It was shown that, under variation 
of the system parameters, abrupt changes of the system qualitative behaviour 
are observed when an impact occurs at a point where the forcing velocity is 
discontinuous (a corner-impact bifurcation point). 

The observed behaviour was explained in terms of appropriate local maps. 
In particular, by using the technique of discontinuity-mappings recently pro- 
posed in [17] and [|12J, it was suggested that a corner-impact bifurcation of 
the oscillator corresponds to a border-collision of a fixed point of the associ- 
ated Poincare map. An important difference was highlighted between corner- 
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Fig. 1: Valve train configuration. 



impact bifurcations and other types of discontinuity-induced bifurcations H) in 
impacting systems such as grazing of limit cycles Il28ll38ll22l|26ll33ll24l . While 
the normal form map of a grazing bifurcation is typically characterised by a 
square root singularity [28J, the local normal form map associated to a corner- 
impact bifurcation was shown to be a piecewise linear map with a gap such as 
those studied in ||T9l . Hence, as explained in ||8l, an appropriate classification 
method needs to be used to investigate this novel class of bifurcations. 

In (T4|. it was conjectured for the first time that corner-impact bifurcations 
are fundamental in organizing the complex behaviour observed in cam-follower 
systems. It was shown that, as the cam rotational speed increases, these sys- 
tems can exhibit sudden transitions from periodic solutions to chaos. Such 
transitions were conjectured to be due to corner-impact bifurcations. 

In this paper, we present a careful analysis of corner-impact bifurcations in 
cam-follower systems. We derive analytically the normal form map associated 
to such a bifurcation in a representative example of interest where the cam pro- 
file is characterised by a discontinuous acceleration. In particular, we investi- 
gate the bifurcation behaviour exhibited by this system under variations of the 
cam rotational speed. We find that following the detachment of the follower 
from the cam, the system can exhibit complex nonlinear phenomena involving 
chattering, period adding cascades and the sudden transition from periodic at- 
tractors to chaos. We explain the sudden transition to chaos observed in the 
system in terms of a corner-impact bifurcation. Namely, we show that dra- 
matic changes in the system behaviour are observed when, under parameter 
variation, one of the impacts characterizing the system trajectory crosses one 
of the manifolds in phase space where the cam acceleration is discontinuous. 

We prove that the normal form map of the corner-impact bifurcation in 
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these systems is a piecewise linear continuous map rather than discontinuous 
because of the higher degree of discontinuity of the forcing signal provided 
by the cam with respect to that of the forcing considered in |8|. We wish to 
emphasize that such a finding is generic for the the wide class of impacting 
systems characterised by forcing terms with discontinuous acceleration. 

As shown in the paper, the derivation of the mapping has an immediate 
practical relevance. In fact, the derivation of a piecewise linear normal form 
map implies that the strategy to classify border-collisions in piecewise linear 
continuous maps due to Feigin [131 can be used, under some circumstances, to 
classify corner-impact bifurcations in continuous-time impacting flows. 

The rest of the paper is outlined as follows. In Sec. 2, we present the mod- 
elling of the cam follower system of our concern where the cam profile has been 
assumed to be characterized by a discontinuous acceleration. Then in Sec. 3 
the numerical bifurcation analysis is presented under variation of the cam rota- 
tional speed. In Sec. 4 we present the analysis of the corner impact bifurcation 
phenomenon detected in the system and we classify the ensuing dynamics by 
using an appropriately derived local mapping. Finally, conclusions are drawn 
in Sec. 5. 



2 Modelling 

The formulation of an appropriate model for a cam-follower system can be a 
challenging task for most applications. Various models with different degrees 
of complexity have been proposed and extensively studied. They range from 
simple models with one degree-of-freedom (DOF) such as that described in 
|20l to complex models characterised by many DOFs, as for example the 21 
DOFs model studied in [36J where additional effects of camshaft torsion and 
bending, backlash, squeeze of lubricant in bearings are included. Nevertheless, 
there is a general agreement in the literature, confirmed by experiments, that 
a lumped parameter single degree-of-freedom model is adequate to represent 
the main qualitative features of the d5mamic behavior of the system of interest 
|3l|20llll[l5l. 

The schematic diagram of the cam-follower system under investigation is 
shown in Figure El We consider the following second order equation to model 
the free body dynamics of the follower away from the cam, 

mq"{t)+hq'{t)+kq{t) = -mg 

ifq{t)>c{t), (1) 

where m, b, k and g are constant positive parameters representing the follower 
mass, viscous damping, spring stiffness and the gravitational constant respec- 
tively. The state of the follower is given by the position q{t) and the velocity 
q'{t). The cam position is given by c ( f ) and we assume that the follower motion 
is constrained to the phase-space region where q{t) > c{t). 

The dynamic behavior when impacts occurs (i.e. q{t) = c{t) ) is modelled 
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via a Newton restitution law as 

q'{t+) = {l+ry{t)-rq'{t-) 

ifq{t)=c{t), (2) 

where (]' {t^) and cj'{t^) are the post- and pre-impact velocities respectively, 
c'{t) is the projection of the cam velocity vector at the contact point along the 
direction of the free movement of the follower, and r G [0, 1] is the coefficient 
of restitution used to model from plastic to elastic impacts. 




Fig. 2: Cam-Follower schematics, (a) t=0. (b) t=T. 

An essential ingredient of the model is the choice of the cam profile, c(t). 
The cam is assumed to be rotating at a constant angular velocity w and can be 
interpreted as the "control action" acting on the follower state as suggested in 
l30l . Therefore, c{t) is carefully selected in applications as a trade off between 
several optimality criteria dependent upon the specific device being consid- 
ered and the unavoidable physical constraints present on the system. 

Typically, this results in a design process where the cam profile is selected 
by using splines and can contain several degrees of discontinuity For exam- 
ple, the cam for a single overhead camshaft valve train is designed by using 
quadratic splines and, as a consequence, discontinuities are present in its accel- 
eration. In general, it is not uncommon in applications, to find cam geometries 
characterised by continuous cam positions and velocities but a discontinuous 
second-derivative (30|. 

In what follows, we assume the cam profile to be characterised by a dis- 
continuous second derivative as shown in Figure |3] For the sake of brevity. 
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Fig. 3: (a) Cam profile, (b) Constraint position c{t), velocity c'{t) and accelera- 
tion c"(i). 



the analytical expressions of the cam profile and its derivatives are reported 
in Appendix A. The case of a smooth cam profile with continuous first and 
second-order derivatives is also of interest in applications and was studied ex- 
perimentally in 0. 
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3 Numerical bifurcation analysis 

The model represented by equations ^ and (|2]l was found to exhibit an in- 
tricate bifurcation behaviour including the sudden transition to chaos under 
variation of the cam rotational speed, cv [14]. The presence of bifurcations and 
chaos was also confirmed by experiments as described in 0. 

Here we briefly summarize some of the most striking behaviour exhibited 
by the system focusing on the abrupt transition from a one-periodic impacting 
solution to chaos observed when w w 673.234445 rpm. 

In general, starting from low values of cv the system exhibits solutions char- 
acterised by permanent contact between the cam and the follower. As cv in- 
creases the follower is observed to detach from the cam during its evolution 
and then to impact with it. A typical periodic evolution with impacts is shown 
in Fig. |4|a) when a; = 183 rpm. We observe that the follower and the cam 
are in contact with zero relative velocity for part of the orbit {sticking) and then 
detach giving rise to impacting behaviour. As shown in Fig|4|b)-(c) a careful 
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Fig. 4: Time simulation at cv = 183 rpm. (a) Follower position, q{t) (Light); 
Cam position, c{t) (Dark), (b) Relative position, q{t) — c(f). (c) Phase 
space, q{t) - c{t) Vs. q'{t) - c'{t). 
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look at the follower evolution shows that a chattering sequence is present, where 
theoretically an infinite number of impacts accumulate in finite time. (Note 
that in practice chattering is associated to a large but finite number of impacts.) 

Chattering can be associated to an intricate bifurcation structure. In Fig.lSJa), 
the location of the impacts in the cam surface is depicted for each value of cu, 
characterising the follower asymptotic solution. We see that following detach- 
ment at about 114 rpm, the follower immediately exhibits multi-impacting be- 
haviour and chattering (characterised by the accumulation of the impact lines 
in the diagram onto the darker areas corresponding to the chattering accumu- 
lation points). An interesting phenomenon is the appearance of resonant peaks 
associated to impact lines crossing the boundaries where the cam acceleration 
profile is discontinuous (represented by dotted lines in the figure). A detailed 
analysis of this bifurcation scenario is presented in {T7\. 

This phenomenon can be classified as due to a corner-impact bifurcation, a 
t5^e of discontinuity-induced bifurcation recently described in H). Namely, 
at certain values of cu, one of the impacts characterising the follower motion 
occurs at a point on the cam profile where the acceleration is discontinuous. 
We shall seek to investigate analytically this phenomenon and classify the be- 
haviour following the corner-impact event in the cam-follower system of in- 
terest. For the sake of simplicity, we focus on a different region of the system 
bifurcation diagram depicted in Fig.|5|c). Here a one-periodic solution charac- 
terized by one impact per period exhibits sudden transitions to chaos as cu is 
decreased below 673.234445 rpm. A close look at the impact bifurcation dia- 
gram in Fig. |5lc) and in the stroboscopic bifurcation diagram Fig. |5|d) shows 
that such transitions occur precisely when the impact characterising the solu- 
tion crosses the cam discontinuity boundaries (the dotted lines in figure |5|c)). 
Specifically, the sudden transition to chaos is due to the corner-impact bifur- 
cation of the periodic solution depicted in Fig. ISje). Past the corner-impact 
bifurcation point the system exhibits chaotic behaviour (see for example the 
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trajectory reported in Fig. |5tf) for cu w 670 rpm). The rest of this paper is 
devoted to the analysis of this bifurcation scenario. 
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Fig. 5: (a) Impact bifurcation diagram for [115,200] rpm. The phase of an im- 
pact (pi (rad), is plotted agaist w. (b) Time evolution for 175 rpm. (c) 
Impact Bifurcation diagram for cv = [660, 750]rpni.(d) Stroboscopic Bi- 
furcation diagram for cv = [660, 750]rpm. (e) Bifurcating orbit at the cor- 
ner impact point at cu = 700rpm. (f) Chaotic evolution for cu = 670rpm. 
Dotted and dashed lines in the diagrams, represent phases where the 
cam profile is discontinuous. Vertical curves in panels (a),(c) shows the 
cam position velocity and acceleration as function of the phase. 
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4 Corner Impact Bifurcation Analysis 

The numerical observations reported above indicate that a corner-impact bifur- 
cation is causing the transition to chaos observed in the cam-follower system. 
Specifically, we are interested in analyzing the occurrence of the corner impact 
bifurcation depicted in Fig. |5jc) when 

w w 673.234445 rpm. Numerically, we detected that the bifurcating orbit, 
shown in Fig. |5|e) is a one-periodic orbit characterised by one impact per pe- 
riod. As the rotational speed of the cam is decreased, at the bifurcation point, 
the impact is observed to cross the point on the cam surface where the cam 
acceleration is discontinuous. To investigate this novel type of discontinuity- 
induced bifurcation we will construct analytically the Foincare map of the sys- 
tem close to the bifurcation point. We will then study the bifurcations of the 
fixed point corresponding to the periodic solution of interest. A crucial point 
in the analysis is to assess whether the resulting map is piecewise linear con- 
tinuous or not. Indeed, only if this is the case, the theory of border-collision 
bifurcations (see l35l[T3ll ) can be used to classify the possible solutions branch- 
ing from the corner-impact bifurcation point |25| . 

We use the concept of discontinuity mapping (or normal form map) re- 
cently introduced in (T7|. | T2| to construct analytically the Foincare map asso- 
ciated to the bifurcating orbit of interest. We use the cam-follower system de- 
scribed in Sec|2]as a representative example to carry out the analytical deriva- 
tions. 

4.1 Poincare Map Derivation 

We are interested in the analysis of the period one orbit at the comer-impact 
bifurcation point. Such orbit is sketched in figure |6l Then, close to such peri- 
odic orbit we define the stroboscopic map P as the mapping from the follower 
state xi e Hi at a stroboscopic time instant ii to the next stroboscopic point 
X2 G T\2. Without loss of generality, we assume that t„ = —^ + {n — l)T 
for n = 1, 2, 3, . . ., where T is the period of the cam forcing cycle (note that 
T = InjuS). Namely, we have: 

X2 = P(xi). (3) 

To construct P we would need to flow forward using the system evolution 
from x\ to X2 for time T taking into account the possible occurrence of impacts 
and therefore applying Newton's restitution law as required. Alternatively, as 
shown in [17J, it is possible to construct P as the composition of three submap- 
pings: (i) an affine transformation Pj jji from the stroboscopic plane TTi at 

t\ = — -J to the plane TTd going through the corner impact point at f = 0; (ii) an 
appropriate zero-time discontinuity mapping(ZDM) Pq on TTo accounting for 
the presence of the discontinuity; and again (iii) an affine transformation P-ij: 11 

from the plane Wd at f = back to the stroboscopic plane 172 at f2 = 7- Specifi- 
cally, while P]^ J j2 and P2 7/2 are fixed time maps that accounts for the follower 
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evolution away from the cam as if no impact had occurred, the ZDM repre- 
sents the correction that needs to be made to the system trajectories because of 
the presence of impacts. Figure |6] represents the global map composition. This 
means that we can write 



P = Pij/i °Pd° Pi,T/2' 



(4) 



where Pj jij '■ Hi h^ T\]j, will map the state from the initial condition Xi on 
the stroboscopic plane TTi to a point xj on the discontinuity plane TTd as if no 
impacts had occurred. Pjj : T\q ^ Ujj will then map x7 to the point xt appro- 
priately correcting the evolution for the presence of impacts (See Fig. [7}. Finally 

P2J/2 ■ ^D 1-^ ^2, 

will map xj" to a point X2 back onto the stroboscopic plane 172. Iri so doing, 
as discussed in |[T7|. (T2|, the effect of the system discontinuities due to impacts 
are all taken into account by the ZDM, Pq, which is therefore often termed as 
the local normal form map in the context of the theory of discontinuity-induced 
bifurcations I 




-T/2 T/2 

Fig. 6: Global map composition. 



4.1.1 Derivation of Pi 7/2 and Pij/i 

As explained above, the maps Pi^r/2 arid P2J/2 ^re defined only in terms of the 
free body dynamics of the follower and the cam rotating period T (depending 
upon the cam rotational speed cu). Therefore we can solve equations ^ to get 
an analytical expression of the flows generating the mappings of interest. 
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Specifically, we define 



Xt 



q'(t) 



yt 



c{t) 
c'{t) 



as the state vector for the follower and the cam respectively. 

Then, as explained in Appendix B, the generalized solution of l|l]l is: 



where C 



2m ■ 



Xt 



Wq 



e ^^ {Icos{wst) + Asin{wst)) xq 

(ptXo, 



(5) 



, Ws 



cUq — C^/ I is the identity matrix, (pt^o 
represents the system flo^v for time t starting from the initial condition xq and 



A 






J_ 



Note that, in general, the system flow operator can be expressed as: 



-Ct 



'Pt 



Ws cos{cVst) + C sin{wst) 



sin{cvst) 



-cuq sin(a;sf) a^s cos{wst) — Csin(a'sO 



(6) 



The submapping Pij/2 can then be easily obtained using ||5| as: 



i,T/2[ 



e-^'^l^ {Icos{cvsT/2) + A sm{wsT/2))x 



4>tx. 

2 



(7) 



4.1.2 Derivation of Pq 

As explained in fl2l , the ZDM can be obtained by an appropriate composition 
of backward and forward flows so that the overall time spent following back- 
ward and forward is zero. As explained earlier, the ZDM is the correction that 
maps the point xj e Hd onto the point x'^ G Hd taking into account the pres- 
ence of impacts in the trajectory of interest. In what follows we assume that 
only one impact occurs over one cycle of the periodic orbit of interest as we 
suppose to be sufficiently close to the bifurcating orbit Xj shown in Fig. [6l Fig- 
ure [Zl shows a schematic diagram that describes the construction of the ZDM, 
close to the corner-impact bifurcations. Without loss of generality we assume 
that the origin is placed at the Poincare section Y\d. To derive analytically the 
mapping xt = Pd(x7) we need to perform the following steps: 

1. Starting from xj, we find the time f, at which the impact occurs. Namely, 
ti is obtained by looking at the difference, ((j(t) — c{t)), between the fol- 
lower position and the cam position close to f = 0. Given a vector z, we 
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indicate by [z]^ its first component. Then q{t) = [xt]^ — ^ and therefore, 
close to xj, tj can be obtained as the nearest solution of the equation: 

H{xZt., ti) := \xZt. - y^t;] ^=h- [(p-tiX^ - y-tj = 0, (8) 

where h = [ 1 ] . 

Hence, tj is implicitly defined by the equation H(xlj , tj) = 0. Once, f/ is 

foxind, the pre-impact state of the system, xlj,, can also be obtained as 

xif. = <P-tiX';^- (9) 

Note that f , can be either negative or positive according to whether the 
impact occurs to the left or to the right of f = 0. 

2. Using the restitution law l|2ll, we can then write the post-impact state of 

r x^j, as 

xtf. = xZt. + R{xZt. - y-t,) = p{xZt., y~t,), (10) 



the follower x^j, as 



where 

[ -(1 + r) 

3. Finally, to obtain xj", we flow forward for time tj starting from the post- 
impact state xitf found at the previous step. In so doing, the state of the 
follower x'^ G T\u can be computed as: 

x+ = cPt,xtj. (11) 

Using equations |[9ll,lfT0ll and lITTI l we can then write explicitly the ZDM as: 

4= Poi^d) =il + 4>t,Rcp~t,)xj-<Pt,Ry-t,, (12) 

with tj defined implicitly by the equation l|8ll. 

4.1.3 Constructing the Stroboscopic Map 

Composing the submappings P-[j/2> Pi.T/l ^^^ ^D given by l[7| and |(T2)| , we 
can then construct the stroboscopic Poincare map, P, of the system close to the 
corner-impact bifurcation point from a generic x„ E T\„ to x„+i e n„+i as: 

X„+i = P{x„,T) = P2,T/2{PD{Pl,T/2{Xn))) 

= 4)T({i + 4>tiR4>-t,)4>TXn-4>t,Ry-t), (13) 



where f , is implicitly defined by the equation H{x„, tj) =h ■ f t _f ,x„ — y-tj 
0. 
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Fig. 7: ZDM construction. 



Note that the fixed point {x* associated to the periodic solution existing for 
a fixed value of the cam period T = T*), can be obtained by solving equation 
l(T3t for x„+i = Xn = X* i.e.. 



I — (pj* + Cpjt^RcpT^ 

2 2 



cpiRyo, 

2 



(14) 



with t* = 0. 

In what follows we are interested in studying such mapping locally to the 
corner-impact bifurcation point detected when cv = cu* = 673.234445 rpm, 
corresponding to a period T* = 0.08912199969159 s. The fixed point associ- 
ated to the bifurcating orbit is x* = [ 5.09700788184250 ] ' . These values 
were detected firstly numerically and then obtained analytically by solving (14)1 
through an algebraic manipulation software (For the sake of brevity we leave 
out the computer algebra here). 



4.2 A locally piecewise-linear continuous map 

Let 6x„ and 5T be sufficiently small variations of the state and parameter from 
the bifurcation point x*, T*. We can then linearize the map x„^i = P{xn, T) in 
l fT3l l about this point as: 



5x„,,='-X^sx„+'-m:nsT. 



dXn 



dT 



(15) 



dP 



For the computation of ^ it is essential to take into account the implicit de- 
pendance of f , on x„ and T. Hence, using implicit differentiation, we have 

dPix„,T) _ dPiXn) , dP{ti)dtiix„) 



dx,i 



dx„ 



dti 



dx„ 



(16) 
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Using ((T3] |, we can then write 

dP{Xn) 



dx„ 



4>t{i + 4>~uR4>u)4>t (17) 



^^ = 4)r (4>'t,R'P-u - {4>tiR4''-t) 4>T_Xn - 4>'t,Ry-t, + 4>uRy'-t) (18) 
Moreover, using implicit differentiation theorem, from (|8]l we have: 

dH{Xn,ti{Xn)) _ dHjXn) ^ dHjtj) dti{x„) _^ 
dXn dx„ dt; dx„ 

The above expression can be used to compute the remaining term in |[16] | as: 

dU{x„) _ fdH{U)\-'dH{x„) ^^^^ 



dXn V ^^i J ^^n 

where 

dH{t,) 



dH{Xn) 



-h-(^4>'r_,Xn-yQ 



-t. ' 



h-4) 



T 



dx„ 2 *!'' 

and /z = [ 1 ] . 

After substituting lfT7b,|[T8]l and iW[ in i fTSt we obtain 

dPiXn, T) 



x„ =x ■ 
T=T* 



<pT* 

1 



dx„ 

( . \ 

I + R) + ({Rcp'o - (P'oR) (Pr*x* + (P'oRyo - Ry'o 



y h ■ U'^.x* - y[, j J 



7 



(20) 

In an analogous way, for the computation of ^, it is essential to take into 
account the implicit dependance of f; on x„ and T. Hence, by using implicit 
differentiation, we have 

dP{xn,T) ^ dP{T)_ dPiU^dUiTl 

Using ( [T3] l, we can then write 
dPjT) _ 



4>'t + \<P\^,R^.,^ + \^^,M'r^ Xn - \'i>\^,Ry-U - '^I+,^%^- 



(22) 
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Again, from @ we have: 

dH{xn,ti{x„)) _ dHiT) , dH{ti)dti{T) 



0, 



dT dT dti dT 

that can be used to compute the remaining term in l(2T) l. Namely, we obtain: 

dt,{T) _ fdH{ti)Y'^dH{T) 



where 



dT 



dH{T) 
dT 



dti 



dT ' 



(23) 



}i- [4>'T_fX„-yLti 



1 i"^A.' ^y-tuT 



and 



dyt,T 
dT 



^c'{t) 



-ic'(0-^c"(0 



Finally, substituting |[T8|, l|221l and ll231l into ||2T]|, yields 

dP{Xn, T) 



dT 



T=T* 



Zl'JZTl/ 2. 2 2 61 



+ cPt* ({Rcp'o - (l^oR) ^i*^* + ^oRyo - Ry'o 



n ■ I ■2(pjtX ^g^ 



h ■ [ (p'j.x* - y[, 



(24) 



We can then compute explicitly these quantities for the cam-follower sys- 
tem of interest. In particular, after some algebraic manipulation, we have: 



A:=ll(.'.T)- 


2 


—r 

L Id-' 


Q+<fd) 

1 




" 
—r 


CpT* 

2 


(25) 


and 




B:-^^(x*,n40.* 


-< + (l + r)c^ 


2 2 






1 

+ 2 


4>T*^ f- 




-4 


^h*d . 


■ (26) 
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Note that both the matrices A and B as defined by (|25)-(|26) depend on the 
value of the second derivative of the cam acceleration Cq at the impact point. 
Therefore the map is actually piecewise-linear locally to the bifurcation point 
where the cam acceleration is discontinuous, i.e. 



c"- := lim c"(t) i^ lim c" it) := c"+. 

Then, the local map can be expressed as: 

A"^x„ + B"<5r, If C-i5i:„ + D-5r< 0, 
A^&Xn + B+iST, If C ■ bxn + D-6T>0, 



6x 



n+l 



(27) 



where 



dx 



B^ 






with the index ± indicating whether the matrices are evaluated with Cg = Cq 
or c^' = c^. 

We have established that close to the corner-impact bifurcation point, the 
d5rnamics of the follower can be studied by means of the local mapping l(27|l. 

Now, from flSl l, the global Poincare map is known to be a continuous func- 
tion of the cam position and velocity through the term i/-t,. Moreover, the map 
is independent from the cam acceleration. It follows, that the map is continu- 
ous at the bifurcation point, i.e. we must have that 



when 

Therefore we have 



A~5xn + B^ST = A+5xn + B+6T, 



C5x„ + D5T = 0. 



C = h-{A' 



A- 



and D = h- (B' 



B^. 



Substituting the numerical values of the map parameters for the cam fol- 
lower system of interest, we obtain the following analytical estimates of the 
map matrices: 



0.82093496821478 
2.52012201452530 



0.01346530915655 
0.82093496821478 



-51.62757990297 
-5455.79455977621 



and 
C = 



0.68571072072040 
2.30988433707948 



-0.07351052377964 
0.68571072072040 



, B^ 



208.11740649865 
-5051.96030903248 



-0.13522424749438 -0.08697583293619], = 259.7449864016200. 
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4.2.1 Numerical Validation 

We will now validate our numerical findings by comparing the map l|27t , which 
was derived analytically, with the numerical estimates of the mapping ob- 
tained by means of simulation and an optimized fitting algorithm close to the 
bifurcation point. 

To derive such an estimate, we use an accurate event-driven numerical al- 
gorithm to simulate the cam dynamics over one period starting from a set of M 
different initial conditions and parameter values. Namely, say 8x„ the vector of 
M possible perturbations of x* and 8T the vector of M possible perturbations 
of T. We then simulate the cam dynamics from each of the perturbed initial 
conditions and parameter values to obtain the vector <5x „+i = x* — x„+i after 
one period. We repeat the set of simulation twice, once with the cam accelera- 
tion set to Cq^ and once with the acceleration set to Cq^. In so doing, we obtain 
numerically the vectors 



^^n+l 



^^l+l 



&X 



n+1 



hr^ 
"^n+1 



We then use a least-squares fitting algorithm to estimate the matrices A- 
and B^ that minimize the error 






{A^ 



5x„ 
6T 



The estimated map matrices found using this numerical strategy are 



A- 



0.82093497830369 0.01346530945739 
2.52012201542191 0.82093496286678 



-51.62757113994 
-5455.79411324739 



0.68571065978423 
2.30988432418263 



-0.07351053029558 
0.68571073479454 



,B^ 



208.11731732063 
-5051.95951604729 



We notice that these numerical estimates are almost identical (up to at least 
5 decimal places) to those obtained analytically earlier in the paper. This vaU- 
dates our analysis and shows the reliability of the analytical derivation used to 
get a leading order estimate of the Poincare map close to the bifurcation point 
under investigation. 

4.3 Classification of the Non-Smooth Bifurcation Scenario 

We can now use the locally derived map (analytical or numerical) to classify 
and explain the bifurcation scenario due to the corner-impact bifurcation de- 
tected in the cam-follower system of interest. In particular, the map derived 
above is a piecewise linear continuous map. As the cam rotational speed is 
increased, the period T of the forcing provided by the cam varies. Correspond- 
ingly, at the corner-impact bifurcation point {8T = 0), the map fixed point 
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undergoes a border collision. Feigin strategy for border-collision bifurcations 
can then be used to classify the corner-impact bifurcation scenario fl3ll . 

The idea is to start by recasting the map l|27|l into a canonical form following 
the procedure presented in [U. Specifically, 

1. We eliminate the term depending on 5T by considering an appropriate 
change of coordinates. In particular if we say ci and C2 the coefficients of 
C, we choose: 



^1 



5x 



6x'^ 



5x\ 
8x}, 






so that the map becomes 



bx 



n+l 



A~5x„ + B5T, If C ■ 5x„ < 0, 
A+6x„ + B5T, If C • 5x„ > 0, 



where 



K - T^d 

1 ci 

2 ci 



■'1 






1525.26226128059 
-615.02768162765 



with Uj: being the coefficients of A . 

2. Then, using the strategy presented in (H [TTj . we consider the change of 
coordinates x = W^^x where the matrix W is obtained as W = T~0~ 
with 



O" = 



C 

CA- 



,T- 



where d-^ is the linear coefficient of the characteristic polynomial of A 
given by p~{X) = A^ + d^A + d^. Applying such a similarity transfor- 
mation, the map matrices become: 



1.64186993642956 1 
-0.64 



A^ 



1.37142144144080 1 
-0.64 



and 



152.75990 
207, 79599 



C= [ 1 ] 



As explained in ||T3l l4|, we can now classify the type of bifurcation scenario 
observed at the bifurcation point under investigation by computing the map 
eigenvalues on both sides of the boundary. For the case under investigation, 
we have that: (i) the eigenvalues oi A~ are A^ = 1.0052 and A2 = 0.6367; (ii) 
the eigenvalues of A+ are A^j = 0.6857 ± jO.4120. Hence, according to Feigin's 
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Fig. 8: Numerical bifurcation diagram of the local map ll27t with the analyti- 
cally estimated matrices. The border collision when ^T = corresponds 
to the corner-impact bifurcation point at cu w 673.2 rpm. Note that as 
predicted a nonsmooth fold scenario is observed with no fixed point ex- 
isting for 5T < and two coexisting fixed points, one stable, the other 
unstable for 6T > 0. 



classification strategy, since the total number of real eigenvalues greater than 
unity on both sides of the boundary is odd, the bifurcating fixed point will 
undergo a nonsmooth saddle node bifurcation and ceases to exist ITSll . This 
is in perfect agreement with what observed numerically as shown in Fig. |8l 
where the local bifurcation scenario observed in the map is shown. 

Therefore, we can explain the sudden transition to chaos observed in the 
cam-follower system under investigation as due to the occurrence of a corner- 
impact bifurcation. Namely, the corner-impact is associated to a nonsmooth- 
f old scenario causing the disappearance of the stable impacting solution under- 
going the bifurcation. This causes trajectories to leave the local neighborhood 
where they are confined before the bifurcation and converge towards the sta- 
ble coexisting chaotic attractor when w is decreased below the corner-impact 
bifurcation point. 

Hence, we can conclude that corner-impact bifurcations in cam follower 
systems can indeed lead to dramatic changes of the system qualitative behavior 
including sudden transitions from periodic solutions to chaos. 

5 Conclusions 



We have studied a novel type of discontinuity-induced bifurcation in a class 
of mechanical devices widely used in applications: cam-follower systems. Us- 
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ing a representative second-order model of the follower, we have shown that 
its dynamics can undergo several bifurcations including sudden transitions to 
chaos as the cam rotational speed is varied. We analysed in detail the corner- 
impact bifurcation of a one-periodic solution characterised by one impact per 
period. In particular, we observed that the system behaviour undergoes dra- 
matic changes when the impact occurs at a point where the cam profile is dis- 
continuous. Using the concept of discontinuity mappings, we derived analyti- 
cally the Poincare map associated to the bifurcating orbit in the case where the 
cam profile has a discontinuous acceleration. Then, using the classification 
strategy for border-collision bifurcations, we proved that the corner-impact 
causes the fixed point associated to the bifurcating orbit to undergo a nons- 
mooth saddle-node bifurcation. Namely, the fixed point ceases to exist, with 
the trajectories being attracted towards a chaotic invariant set. 
We wish to emphasize that: 



• 



• 



the analysis presented above applies with minor changes to the case of 
impact oscillators forced by signals with discontinuous second deriva- 
tive. As shown above, this leads to maps which are locally piecewise lin- 
ear continuous close to a corner-impact bifurcation point. This extends 
the analysis presented in [8J for the case of an impact oscillator forced 
by a function with discontinuous first derivative. We conjecture that the 
properties of the local mapping depend on the degree of discontinuity of 
the forcing signal. This is the subject of ongoing work. 

As shown in |12j, discontinuity-induced bifurcations in flows are usually 
associated to maps which are not piecewise linear. Grazing bifurcations 
of limit cycles where known to be associated to maps with square-root 
singularities in impacting systems and Filippov systems 1261 or maps 
with higher order nonlinear terms in the case of piecewise-smooth con- 
tinuous flows 
(PWSC). The only cases in the literature where the map was indeed found 
to be piecewise linear-continuous were corner-collisions in PWSC sys- 
tems and grazing sliding bifurcations in Filippov systems. So far, no evi- 
dence was given of a bifurcation event in impacting systems associated to 
locally piecewise-linear continuous maps. The corner-impact bifurcation 
scenario presented in this paper fills this gap in the literature. 

We believe cam follower systems are a particularly useful set-up to show 
generically the behaviour of impacting systems with discontinuous forc- 



Finally, the results presented here can open a way to future work towards a 
better understanding of the complex dynamics of cam-follower systems. This 
can lead to less conservative solutions to detachment avoidance, hopefully 
without recurring to highly stiff closing springs and maybe active control strate- 
gies. 



5 Conclusions 



22 



Acknowledgments 

The authors wish to thank the anonymous reviewers whose comments led to a 
consistent revision of the original version of this manuscript. They also grate- 
fully acknowledge support from the European Union (EU Project SICONOS - 
V Framework Programme, IST2001-37172) and the project MIUR-PRIN MACSI 
funded by the Italian Ministry for Research and University. The paper was 
completed during a research visit of the authors at the Centre de Recerca Matem- 
atica in Barcelona thanks to support from the Government of Catalunya. 

Appendix A - Cam Profile 

We report below the analytical description of the representative cam profile 
considered in this paper. As shown in Fig. |9l in this case the cam profile is the 
result of a geometrical based design. 
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Fig. 9: Cam profile definition, (a) = 0. (b) = n. . 

The lift profile c{0) can be defined from the construction as a piecewise 
smooth function of the angle as: 



c{0) 



coid) 


IfO< 


< f - 01, 




cm 


Iff- 


01 < < f - 


-02, 


C2{0) 


Iff- 


02 < < f - 


-03, 


cm 


Iff- 


03 < < TT, 





(28) 
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co(0) = Po 

ci(0) = -(Cisin(0 + 0i) + (pf-K?cos(0 + 0i)^)3 

C2(0) = K2 sin(0 + 03) + (p^-K^cos(0 + 03)2)5 

C3(0) = P3 

where 0i = ZSOR, 02 = ^SOP, 83 = ZSOQ. Additionally, k, and p, are 
constant parameter given by our particular geometrical construction of the cam 
as {See Fig|9) 



Ki = \\qOi\\, Po = \\qR\\, P2 = IIO2PII, ^29^ 

K2=\\002\\, Pi = ||OiR||, P3 = ||OQ||. ^ ^'' 
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